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SUMMARY 


Linearized theory for corapressible unsteady flow is used to derive 
the velocity potential and lift and moment coefficients in the form of 
power series in terms of the frequency of oscillation for a haarmonically 
oscillating rectangular wing moving at a constant supersonic speed. 

Closed expressions for the velocity potential and lift and moment coef- 
ficients aasociated with pitching and translation are given to the sev- 
enth power of the frequency. These expressions extend the range of use- 
fulness of NACA Report 1028 in which similar expressions were derived 
to the third power of the frequency of oscillation. For example, at a 
Mach number of 10/9 the expansion of the potential to the third power 
is an accurate representation of the potential for values of the reduced 
frequency only vp to about 0 . 08 ; whereas the expansion of the potential 
to the seventh power is an accvtrate representation for values of the 
reduced frequency up to about 0.2, a value of this parameter large enou^ 
to cover most rectangular -wing flutter cases likely to occur at this Mach 
number. 

The section and total lift and moment coefficients are discussed 
with the aid of several figures . In addition, flutter speeds obtained 
in the Mach number range from 10/9 to 10/6 for a rectangular wing of 
aspect ratio ^<-.53 by \islng section coefficients derived on the basis of 
three-dimensional flow are compared with flutter speeds for this wing 
obtained by using coefficients derived on the basis of two-dimensional 
flow. 


INTRODUCTION 


A method for obtaining the air forces and moments acting on harmon- 
ically oscillating rectangular wings in supersonic flow is given in ref- 
erence 1. This method is based on the expansion of the velocity potential 
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in powers of the frequency of oscillation and is applied to harmonically 
translating and pitching rectangular wings to obtain expressions for the 
associated forces and moments in^volving the frequency to the third power. 

In the expansion of the velocity potential, the frequency of oscil- 
lation and Mach number enter the results in a combined form such that 
the range of frequency for which the expansion to a given power repre- 
sents an acceptable approximation to the actual potential decreases as 
the Mach mmiber decreases toward unity. As a result, the force and 
moment expressions of reference 1 apply to a sufficiently broad frequency 
range for most flutter studies only if the stream Mach number is greater 
than about 1.5- In order to obtain results that cover a larger part of 
the transonic range, that is, results that apply over a wider frequency 
range at Mach numbers nearer unity, the present paper extends to the sev- 
enth power of the frequency the expressions of reference 1 for velocity 
potential, section force and moment coefficients, and total force and 
moment coefficients . This extension results in a coverage of frequency 
that is generally sufficient at Mach numbers as low as about 1 . 1 . It 
also, of course, increases the frequency range covered by the approxi- 
mate theory at all supersonic Mach numbers . 

Althou^ the method of reference 1 furnishes a straightforward 
approach for malting this extension, the present paper employs a more 
concise method based on the velocity potential for a semi-infinite wing, 
developed in reference 2 . In this reference Stewartson makes vise of the 
Laplace transformation to obtain the potential in the form of a definite 
integral. The integrand of this integral expanded in powers of the fre- 
qxiency and integrated term by term can be made to yield results that are 
IdenticaJ. in form with those obtained by the method of reference 1 . For 
the sake of completeness Stewartson *s derivation of the velocity poten- 
tial for the semi-infinite wing is reconsidered herein. 

For illustration the extended section coefficients are used to cal- 
culate flutter characteristics for a rectangular wing of aspect ratio 4.55 
at several Mach numbers in the low st5>ersonlc speed range. These results 
are compared with calculations made by using two-dimensional aerodynamic 
coefficients, obtained from reference 5* 

SYMBOLS 


A aspect ratio 

b one-half chord 


c 


speed of sound in undisturbed medium 
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5’n(x,y 

h = F 

h 

• • 

h,a 

Jo((y) 

k 

Li, Mi 

Li^Mi 

M 

Ma 

P 


section lift coefficient 
total lift coefficient 

section pitching-moment coefficient about axis of rota- 
tion Xq 

total pitch i ng-moment coefficient about axis of rota- 
tion Xq 

f motions defined after equation (17) 


n(i^y) 


verticel displacement of axis of rotation Xq, positive 
downward 

an^jlitude of vertical displacement 

time derivatives of h and a, respectively 

Bessel fmctlon of zero order (first kind) 

reduced frequency, oib/v 

con^jonents of section force and moment coefficients, 
respectively, defined in eqmtions (29) to (52 ) } 
i = .1, 2, 5^ and, k 

components of total force and moment coefficients, 
respectively, defined in equations (37) ani (38) J 
i = 1, 2, 3^ and, k 

Mach number, V/c 

aerodynamic section moment on wing about axis of rota- 
tion Xq, positive leading edge up 

total aerodynamic moment on wing about axis of rota- 
tion Xq, positive leading edge up 

local pressure difference 

aerodynamic section normal force, positive downward 
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p 

total aerodynamic force on wing, positive downward 


nondimensional radius of gyration of wing section about 
elastic axis, ^ where is mass moment of 

inertia per unit span about elastic axis and m is 
mass of wing per unit span 

s 

one-half span of wing 

t 

time 

V 

velocity 

w(x',yniS‘t) 

vertical velocity at surface of wing along chordwise 
section y' = y^' 

Xa 

location of center of gravity of wing measured from 
elastic axis (see ref. 3) 

Xo 

abscissa of axis of rotation of wing (elastic axis) 

x,y,z 

nondimensional rectangular coordinates attached to wing 
moving in negative x-direction, referred to wing 
chord 2b 

x’ = 2bx 


II 


z ' = 2bz 


Zb 

first bending mode shape of wing 

Zta 

vertical displacement of any chordwise section of wing 

2b 

first torsion mode shape of wing 

a 

angle of attack, positive leading edge up 

°b 

effective angle of attack due to vertical translation, 

fi/v 

°D 

amplitude of angle of attack a 


p = \ju^ - 1 
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^ = 2y/A 


0h 

phase angle 

between section lift due to h 

and 

velocity 

h 


0a 

phase angle 

between section lift due to a 

and 

position 

a 


®hm 

phase angle 

between section moment due to 

h and 

velocity 



V 

phase angle 

between section moment due to 

a and 

position 

a 


®A 

phase angle 

between total lift due to a and 


position 

a 


®AM 

phase angle 

between total moment due to a 

and 

position 

a 



K density parameter, jtpb^/m 

p density in undisturbed medium 

(f) disturbance-velocity potential 

CD frequency of oscillation 

first bending frequency of wing 

0 ^ first torsion frequency of wing 

H5 = 2M2/p2 


ANALYSIS 

Velocity Potentials for Harmonically Oscillating 
Rectangular Wings 


As a first step in the analysis, an Integral expression is developed 
for the velocity potential for a harmonically oscillating semi-infinite 
rectangular wing (fig. 1) . As mentioned in the introduction this expres- 
sion is given in reference 2 and is redeveloped herein for the sake of 
completeness . Prom the expression for the semi-infinite wing, the poten- 
tials for the various regions (see fig. 2(a)) of a finite rectangular 
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wing are obtained either directly or by appropriate modifications . The 
analysis for the finite wing, like that of reference 1, is restricted 
to the condition that the Mach line from the foremost point of one wing 
tip does not intersect the opposite tip ahead of the trailing edge . 

Velocity potential for semi- infinite rectangular wing .- Consider a 
thin, flat, semi-infinite, rectangular wing moving at a constant super- 
sonic speed in a chordwise direction norma!!, to its leading edge as shown 
in figure 1. The differential eqiiatlon satisfied by the disturbance- 
velocity potential for the wing ('vdien referred to a rectangular coordi- 
nate system x',y',z' moving tmiformly in the negative x' -direction 
with the x'y' -plane coincident with the mean position of the wing) is 


o2yat ix'j ^,2 3y,2 


( 1 ) 


where x* = 2bx, y' = 2by, z' = 2bz, 2b is the wing chord, and c 
is the speed of sovmd in the undisturbed medium. The boundary condition 
of tangential, flow at the surface of the wing, in accordance with small- 
disturbance linearized theory, can be expressed as 



>■±0 


w(x*,yia’,t) 


V 4. 

^t 


( 2 ) 


where ^ is the vertical displacement of the ordinates of the surface 
of any chordwise section of the wing stich as y' = ym' figure 2. 

The wing is assumed to be executing simple harmonic motion with respect 
to time t, so that t enters only in the exponential exp(iojt), where 
(D is the frequency of oscillation. Equation (l) thus becomes 


= q 2 4. 2iuM ^ ^ . 

ax-2 c Bx- c2 


( 5 ) 


where the disturbance-velocity potentieil 0 is related to by 


0(x’,y',z’,t) = ilf(x',y',z')eiiiJt 
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For the case of the semi-infinite wing having identical motion in every 
chordwise section the boundary conditions that equation (3) must satisfy 
are : 


) = V(x') 

^^7z'=±o 

(x' 

,y' ^ 0) 


loulf + V ^ = 0 
Sx' 

(y' 

< 0) 

(5) 

ijr = 0 

(x’ 

< 0 ) 

(6) 

0 

t 

(z’ 

± 00 ) 

(7) 


Equation (4) is derived from equation (2) and implies that the normal- 
velocity distribution on the wing is given; equation (5) is the condi- 
tion that the pressure be zero off the tip of the wing; equation (6) is 
the condition that no disturbances be propagated forward of the wing; 
and equation (7) is a condition on the behavior at Infinity (the manner 
of approaching zero is associated with the radiation condition of 
Sommerfeld) . Equations (5) to (7) constitute the boundary-value problem 
for the velocity potential 0. 

Applying the Laplace transform 




= r 
0 


e“®^ T];(x',y',z')dx' 


to equations (3) to (7) yields the transformed boundary- value problem 
in the form 


2aM 


/ p p ^ cUm tEr\-r _ 2r- 
— V + — V = lP^s2 + i s - Hr = 


(8) 


= w(s) (y* ^ 0) (9) 

/z'=±0 
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t = 0 (y’ < 0) (10) 

Tjf >0 (z' > ±00) (11) 

The problem in this form is, an pointed out in reference 2, similar to 
the problem treated in article 308 of reference By applying the 
method of reference it- to the present problem, it can be shown that the 
value for ijr at the upper surface of the wing (z’ = + 0 ) is given by 



For convenience equation (l2) may be rewritten (see, for exan^tle, p. 478 
of ref. 5) as 



If use is made of a table of Laplace transforms (for example, pairs 55 
and 89 of ref. 6) and the Paltung or convolution theorem, the inverse 
of the Laplace transform in equation (15) may be written as 


where 


X ^ 

- - i r w(x' - |')G(S',y’)a|' 

p Jo 




G(x',y') 



(x‘ < 3y') 


G(x',y') 




- ^S^'^sec\^dTl 


-i — x' 

. cp2 


x' > 3y') 
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and I' = 2b|. From equation the -velocity potential in terms of 

the nondimens ional coordinates x, and y becomes 


0(2bx,2by,t) = _ 2b^f«t J w|^(x - G(2b|,2by)d| (15) 

■where 


G(2bx,2[by) = jQ^g-x^e"^ (x < 3y) 


G(2bx,2by) = 


(x > Py) 



CO = 

vp2 


Equation (I5) is the desired Integral expression for "the "velocity poten- 
tial for a semi-inf ini-te -wing. Note that for x < py equation (15) 
reduces to the potential for a 'two-dimensional wing (see ref. 5) . 

An expansion of the Integrand of eq'uatlon (15) in po'wers of the 
frequency of osci lla tion is en^ployed because the indicated Integration 
does not appear -to be obtainable in terms of kno-wn functions . Q?he result 
of "this expansion is 


0 =- 


l^be^ 

itp 


X ^ 

jT w[2b(x - I)] ^(-icolf 


q=0 


r=0(q - 2r) 




(16) 
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where denotes the integral part of q/2 and 


Ir(x,y) = J 


(x < py) 


Ir(x,y) = I 



^ sec\)'' dTi 


(x > py) 


It may he shown by induction that 3^ for x > Py can be expressed in 
terms of the functions Fn of reference 1 as follows : 


Ir = 


2r 

^2r 




nt=0 


where 


Fn(xjy) 



x^“ls 


dx 


Two relations involving Fj^ that are of particular importance in the 
next section are 



x“F^(x,y)dx 


m + 1 \ 


(18) 


and. 


Fn(x,y) 



1 

n 


^ ajjjX^“”i(Py)“"^ 


n^l 


\fc-(x - Py) 


( 19 ) 


where 


ai - 
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_ 2“"^(n - l)(n - 2) . . . (n - m + 1 ) /_ ^ .n 

am - ■; : liii > 1) 

(2n - l)(2n - 3) ... (2n - 2m + l) 


Equations (16) to (19) provide the means for obtaining the potential to 
any desir*ed power of the frequency of oscillation. They are used in the 
next section to develop expressions to the seventh power of the frequency 
for the velocity potentials for the various regions of a finite rectan- 
gular wing undergoing torsional and vertical translational oscillations. 
These expressions are extensions of similar expressions developed in ref- 
erence 1 to the third power of the frequency. 


Velocity potentials for finite rectangular wing .- The coordinate 
system and regions of interest for the finite rectangular wing are shown 
in figure 2. On the portion of the wing between the Mach cones emanating 
from the foremost point of each tip (region N in fig. 2(a)) no inter- 
action takes place between the flow on the upper and lower surfaces of 
the wing. On the portions of the wing within the tip Mach cones 
(regions Tp, T2, and T^ in fig. 2(a)) interaction does take place between 
the flow on the ipper and lower sxirf aces . The velocity potential at a 


point in one of these regions is designated by or 

according to the region that contains the point. 


For the particular case of the wing independently performing small 
sinusoidal torsional oscillations of amplitude cXq about some spanwlse 

axis Xq and small sin\isoidal vertical translations of amplitude hQ, 
the equation for is (see fig. 2(b)) 


2^ = ei^ |2bao(x - x^) + h^ = 2ba(x - Xq) + h (20) 

Substituting this expression for into equation (2) gives 


w(2bx,t) = Va + 2b&(x - Xq) + h (21) 

Equation (21) conplles with the restriction Inherent in the development 
of equation (15) , namely, that the wing have Identical motion in every 
chordwise section. 

The velocity potential for region Tq for the wing motion represented 

by equation (20) is obtained directly from equation (16) . Prtting the 
expression for w(2bx,t) given by equation (21) into equation (16), 
retaining only teimis in the expansion ip to the seventh power of cii, 
insertiag the appropriate values of I3, obtained from equation (17), 
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performing the necessary integrations with the aid of equation (iP,) 
yields the folOo«lBg form for the velocity potential to reS^ “ ' ' 


=-a (£ 


where 


p |( h + Va)Ai(x,y) + 2ba[A2(x,y) - XqAi(x,:^ 

J 


(22) 


^l(x,y) = 


- imF, 


0^ 


^ ‘ ^ + (232 - l)F^ + 

£3(1j2 + 5)p2 - 3xF4 + 23%^ . ^ _ 

^ 96 om^^ 


30x2f 4 + (set . tp2 + 3)p^ . -^|t5 (sp't + 28p2 + 2 i)f2 . 


5?6om^' 


10x5 (6p2 + yjj. ^ ^ gp6j, 


155^ 


7x' 


8 o 61 iom^L 


(83^ 


3632 + 33)p2 - 1X)5x^(232 + 3^F4 + 105x2f6 + 

(163^ - 83^ + 632 _ 5^F^ ■ 

nX 

^2(^^y) = / Ai(x,y)dx 

The integrations required to obtain the function Ag may be readily 
performed with the aid of equation (18) . The Integrated values of Fn 
(from eq. (19)) needed In equation (22) are listed in appendix A. 

Use is now made of the potential 0^. equation (22), to obtain the 
potentials for regions N, T2, and T3. The potential ^ is obtained 
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from by substituting y = x/g in the fmctions A]_ and 7 ^ in 

equation (22), in -which case Fn becomes 

iT 


and the -velocity potential (f^ is gi-ven by 




(h + Va)fi(x) + 2 bd|f 2 (x) - i ( 25 ) 


■where 


^,(x) = X - f - . 5)f . I^(2p2 . 5 )^ . . 


5 12m2 


^ 192M^ 


l«,j2 35)£ . + 56p2 + 63)=^ - + 

^ 5 96om^ ^ ' D 11520MP ^ 

5 T 8 p 2 + 251) Y + — (163^ + 2l6p^ + + ^29)^ 

f2(^) = [ fl(x)dx 

Note that because of the -two-dimensional character of the problem in 
region N the spanwise -variable y is not contained in Equa- 

tion (25) can also be obtained by expanding to the se-venth power of c 5 
the -velocity potential for the -two-dimensional -wing gi-ven in reference 5* 
The potential obtained from by replacing y in the func- 

tions Aq and 7^2^ in equation ( 22 ) by A - y, where A = 2 s/ 2 b is the 
aspect ratio. The potential in region T^ (this region exists if 

1 < A3 < 2) is a simple s-uperpcsition of the potentials for regions N, 

Tq, and T2^ as indicated in reference 1, and may be written as 

9^5 = ^^q + 0T2 " 



lil- 


NACA TN 5076 


Forces and Moments 

Section forces and moments The expanded velocity potentials for 
the finite wing are now used to obtain expressions for the section forces 
and moments at any spanwise station of the wing. Since the distribution 
over the entire wing is symmetrical with respect to the midspan section^ 
only expressions for the forces and moments at any station of the one- 
half span adjacent to the origin (see fig. 2) need be considered. 

The local pressure difference between the upper and lower svirfaces 
of the wing may be written as 


The section force, positive downward, is therefore 



(25) 


(26) 


and the section moment, positive leading edge vip, about the arbitrary 
axis of rotation x = Xq is 


Mq^ = J' (x - dx (27) 

IMder the previo^lsly mentioned restriction that the Mach l^ne from 
one tip not intersect the opposite tip ahead of the trailing edge, eq.ua- 
tions (26) and (27) must be evaluated for the two cases that can arise 
(see fig. 5) • These canes are: (l) The Mach lines from the tips do 

not intersect on the wing (Ap > 2) and (2) the Mach lines intersect on 
the wing but the Mach line from one tip does not intersect the opposite 
tip ahead of the trailing edge (l Ap 2) . Only the fi n al forms of 
the section forces and moments are given. These forms are calculated by 
deriving the pressure difference for the different regions from the 
appropriate velocity potentials, making use of figure 5 "to determine the 
limits of integration for the regions involved, and using the relation 
given in equation (18) to perform the integrations involving the func- 
tions F^; the results can be written as 
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P = -4pbV^^e^ 


b V" 




Mq, = -4pb^V^k^e^^ 


^Mi + IM2) + oc^^Mj + 


(28) 


The qiiantitles L^ (l = 1, 2, 5, and 4) and Mi (i = 1, 2, 3 , and 4) 


are the components of the section force and section moment coefficients, 
respectively. The reduced frequency k is related to cd and o5 by 
the relations 


k = :^= ^ 53 
V 2 M 2 


In order to obtain self-consistent expressions for the forces and moments, 
the terms of the velocity potentials associated with h and a are 
expanded to the seventh power of 53 and the terms associated with a 
are expanded to the sixth power. As a result, the coefficients in equa- 
tions (28) are as follows: 

Case l(see fig. 3(a)): For any section between the tip and the 

point where the Mach line intersects the trailing edge, or where 
0 < y < l/p. 


Ll 


+ (232 + ijpj - I 3F2 - 63^5 


1 

24o 


5(63^ + 7)^2 - 50 F^ - + 5(163^ 


+ 



1 

20160 



+ 4832 + 


53) f 2 - 


35(83^ + 9)f 2^^ + 105Fg - 1123 + 


/-\6 


^12838 + 6436 - 163^ + 83^ - 5^F0 


(29a) 
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L2 = 


«pk 


+ 


(2p2 _ i)f 2 - 5P^3 


''iLf + -1 . 
,mJ 2k 


5 F 2 - (sp^ - k^^ + 3 )^ 1 ,. + 


lOp^F^ 


'c5\^ 1 


\MJ 2880 
8p^ + 6p2 - 5)55 + 56p^y 


2 l( 2 p 2 + 3)f2 - TOF 4 - 5 (l 6 p^ - 




(29b) 


L 5 = 


It 


PM%si^ 


2(p2 + i]fi - p\ - (6p^ + 5p2 - i)f 2 + p2(6p2 + 5 )f^(0 


i j^p 2 + 5 ) F 2 + 12P^2 - (lK)p^ + 20p^ - 5P^ + 5 )fi^ 


2p^(l5p2 + 13 )f 5 ^ 3(2P^ + 2ip2 + 21) F2 - 

lo(p2 + 7 )Fi^ + I20p%5 - 3(u^P® + 56p^ - lij-p^ + 


7p2 - 5 )f 6 + 8p^(28p2 + 25) F^ 


f) 


(29c) 


Lj^. = 2P^1 - (3p2 + 1)F2 + I 3(2P^ - p2 + i)f 2 - 2k^^F^ + (20p^ + 


7P^ - 1 


(o\^ 1 


m; 240 11' 

24op^^ + (l68p^ + 64p^ - up2 + 3 ) 1 ^ 


5(p2 + 7 )f 2 + 10 (8p^ - 4p^ + 3P^ - 3 )f 4 - 


ISV + 1 

Im; 20160 


7(2p^ + 27p2 + 


33 )f 2 - 35 (p^ + 9 )fi^. + 2i(i6p® - 8p^ + 6p^ - 5p2 + 5)Fg - 896p^Y + 

(29d) 


( 576 P® + 232p6 - 44p^ + 15P^ " 


f-f f - 
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Mn = ^ 






+ (3p2 + 2)f3 - I 2F2 - 3P^F3 + P^( 5 P^ + i^)F5 


'isf + 


2^ 

30 


!(6p2 + 7) F2 - 5Fi^ - 5P^5 + P^(7P^ + 6 )f^ - 


504 o 


3(l6p^ + 48p2 + 33) F2 - l4(8p2 + + 35Fg - 


a8^^y + ^S(jp2 + - 2x„Li (30a) 

“2 ' [®2 - - ^)*2 + + K§5f2 - ^(8^'‘ - 

itp2 + 3)1^ + (altplt + 8 j 2 . l)pg](|-f - ^|i5(2^2 + 3)|g . 10^^ . 


M/ 2880 


3(l6p^ - 8p^ + 6p2 - 5 )f 5 + (6ltp^ + 24p^ - l^p2 + - - 2 xqL2 


(30b) 


Mt; = 


1 


jtpM^^ 


j^^pS + l^Fg - I ^pS-i - 3(6p^ + 3P^ - 1)^2 + (20P^ + 


21p2 + 3 )fi^ i 5^p2 ^ ^ 32P%5 - 2(l|Op6 + 20P^ - 


5p2 + 3 )\ + (56p6 + + up2 _ l)Fg (l) - + 

2ip2 + 2 i)f 2 - 15(p^ + t)fi^ + l6op^^ - 3(u2p® + 56p^ - 

(Equation continued on next page) 
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+ 73^ - 5)^6 + (1923® + 2323^ + 1 ) 43 ^ - 532 + 


2Xq^Mi + + 2X0L2) 


(50c) 


M4 


= — ^ |_8|^3^3^ + (232 + i)f^ + ^ (23^ - 32 + i)f2 - 63^5 + 
32(53^ + + t)f2 + 5(83^ - ^3^ + 53^ - 5)f]^ 


903^5 + 23^(2832 + 19) Fy 3(23^ + 2732 + 33) ?2 - 


50ll0 


14(32 + 9)fi^^ + 7(163® - 83® + 63^ - 532 + 5)ig - 22l^3^Y + 


83® (153^ + 11)^9 - 2xo(M2 + L 4 + 2X0L2) 


(30d) 


■wiiere Fjj^ (n = 1, 2, . . . 9) ^ given in appendix A, is F^ evaluated 

at X = 1. For any section between the point where the Mach line inter- 
sects the trailing edge and the midspan, or where I/3 < y ^ A/2, 


Ll 



+ 5/d5f ^ 
2it 



643® + ^323^ + 7923^ + 429/ 05^^ 
322560 \mJ 


(31a) 


L2 


2_ 1 + ^3^ + T / caf ' 83^ + 3632 + 35 /a 5 \® 

3k 5 \m/ 192 \ kJ 11520 Vm/ 


(31h) 
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Lx = 


|b2 + 1 _ + 33B^ + 33/ c5\^ 

pM 2 fc 2 [_ 12 Im/ 960 '.M; 


83 ^ + l40p^ + 3373 ^ + 231| 

8 o 64 o 




- 2xoLl 


L4 = 


3^k 


r 2 _ 1 + + 5[ (a\^ 43 ^ + k9B‘^ + 6^ /aS\^ 

2k \uj 2880 W 


83^ + 1803^ + 5613^ + k29( aS\^ 
322560 Vm/ 


- 2 xqL 2 


Ml = — 


1 

B5 


4 _ 43 ^ + ^ 83 ^ + 283 ^ + 2 1/53 

3 15 W/ 560 




643^ + 4323 ^ + 7923 ^ + 4 29 

l 8 l 44 o 


1 / 


- 2 xqLi 


^2 = ^ 
3k 


1 _ 3 tef ^ 3(43^ + 7 )fe \^ _ 7(83^ -t- 3632 + 33) 
576 \M/ ~ 46 o 8 o 


'-\6 


% = 

3M"^k2 


j2 + 1 3^ + 3 M^ ^ 43 ^ + 333^ + 33/ fe\^ 

8 W i76 


83^ + i 403^ + 3373^ + 23 l fm'^ 
46080 W/ 


(31c) 


(31d) 


(32a) 

- 2xoL2 

(32b) 


- 2 xq(Mi + L 5 + 


2xqLi) 


(32c) 
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M,. = 


1 

4(p2 - 1) 

p5fc 

5 


05 


+ 63/^co 




15 


8p^ + 180P^ + 56 ip^ + 429 /55' 
I81IAO 


1680 


t) 

— 2Xq^M2 + + 2xqL^ 


(52d) 


Caise 2 (see fig. 5(b)) : For any section between the tip y' = 0 

and the point where the Mach line from the tip at y' = 2s intersects 

the trailing edge ^r where 0 < y < A - the components of the sec- 
tion force and moment coefficients are given by equations (29) and (30 ) , 
respectively. For any section between the point where the Mach line 
from the tip at y' = 2s intersects the trailing edge and the midspan 


or where 


A-±<y< 


A/2) 


the components of the section force and 


moment coefficients, respectively, are obtained by first adding to equa- 
tions (29) and (30) the results of siibstituting A - y for y in these 
equations and then substracting equations (31) and (32) from the result. 
For example, the term of Lj_ that does not contain o5 is given by 


JL 


^ |-4p2[Fi(y) + F^(A - y]j + ^2p2 + l) " y]| ^ - 

(55) 

Total forces and moments .- Expressions for the total forces and 
moments can be derived by considering o nly case 1 of the previous sec- 
tion. It can be shown that case 2, although more cumbersome to handle, 
leads to the same expressions. Therefore the total force, positive 
downward, may be written as 


pl/p pA/2 

P = 4b / (P)^^ dy + 4b J (P)jj dy 


0 


(54) 


and the total moment, positive leading edge up, about the axis x = Xq 
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In eq.uations ( 3 ^) and (35) "tlie quantities 


(P)^^ and are the 


section force and moment, respectively, whose components are given in 
equations (29) and (30), and (P)jj and (Mo,)n are the section force 

and moment, respectively, whose con5>onents are given in equations (31) 
and (32) . 


Upon performing the Integrations indicated in equations (3^) and 
(35) with the aid of the relation 


r ? ay = S(2_i_S). JL 

Q n(n + 1) 83 


the results can he written as 


where 


P = -8ph%%2^^ 

^ = -Sph^v^k^Ae^ 






L-, = -i. « 1 _ + QP^ 28p^ + 21/ co’) _ 

o3 2h [uj 9^ \M/ 


(36) 


64b^ + Jj-32B^ + 792S^ + k29l (j^^ _ 1 


322560 


W 3A^ l!^ 


2 + 2 P^ + 8P^ + + 


20 


.m; 


+ l8p^ + + 32 M^ p 8 ^ ^ 2 p 6 ^ 2603^ + 256p2 + 128 ^\^ 

W " 6o48o W/J 


840 


1 


(37a) 
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2Xo(mi + L5 + 2XoLjj 


(58c) 
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M4 


1 . - 1) ^ _ li-B^ + 49B^ + 65 /s\^ ^ 

p5k I 5 15 \My 1680 [mJ 

8b^ + i8ob^ + ■ 56 ib^ + 429 ^ JL 

l8l¥w yyy 3 AB 1_ 3 

3 P^ + i6b2 + i6/ffi\^ + lop^ + 2kp,^ + 


2ko 


M 


Ij-aDO 


|3 _ + 

1 - 


2xq(^M2 + L4 + 2XQL2) 


(38d) 


SOME PROPERTIES OF FORCES AMD MOMENTS 


Examination of some of the properties of the extended section and 
total forces and moments of the present paper may be of Interest. For 
this purpose consider, first, the conventional coefficients for section 
lift c^ and section pitching moment c^ that may he obtained from 

eqmtions (28) as 


C2 = - 


= k 




phV2 

_ Mg ^ _2 

2pb2v2 [_ 


-ik(Li + 11.2) °h + ^^(^3 + ^k)°- 
-ik(Mi + iM2^aj^ + + iM^^ja 


( 59 ) 


where = h/v is the angle of attack due to vertical translation. 

From equations ( 39 ) the section lift- and moment-curve slopes (complex 
derivatives) associated with vertical translation and pitching are, 
respectively. 



^ = i2k(Mi + IM2) (iK)a) 


-iI*-k(Li + 1X2) 
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and 


^ + 114) 


^ = -a^(M5 + 1M4.) (tob) 


The coefficients c^ and Cj^ and their associated slopes reduce to 

well-known results in the steady case (k = O) . For exan5)le , substitu- 
tion of the components given in equations (31) into equations (39) yields 
the result C2 = 4a/3 (Ackeret's result) for a nonoscillating two- 

dimensional wing. Eqmtions (39) and (40) have the additional feature 
that multiplying the components and (i = 1 and 2) associated 

with vertical translation by k and the conqponents and Mj^ 

(1 = 3 and i^) associated with pitching by k^ makes these quantities 
more uniform in magnitude and therefore more suitable for plotting. 


Now consider the coefficients for total lift a“d total moment 

Cjn that may be obtained from eqmtions (36) as 


Oi <= - = 41^1^(10. + + k2(L3 + 114)^ 

Cn = ^ - -2[Ilk(Mi + ji^)oh + k®(Mj + IMk)^ 


(la) 


From eqmtions (i|-l) the total lift- and moment-curve slopes (con^jlex 
derivatives) associated with vertical translation and pitching are, 
respectively. 


and 

cLCl 

da 


-i4-k(^Li + il.2) 


kk^^Lj + llj 


dc% 


i2k^M3^ + 



(42a) 


^ = -2k2^M5 + iMi^) (42b) 


Eqmtions (ItO) may be med to illustrate the extent to which the 
approximate expressions of the present paper for the section force and 
moment on a finite rectang\ilar wing may be meful. A comparison, based 
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on these equations, is made in figure of some exact and approximate 
calculations for a two-dimensional wing vertically translating and 
pitching about its midchord position. Results for a two-dimensional 
wing are chosen because in this case approximate and exact results can 
be compared and the extent of convergence determined . The two-dimensional 
caee is a legitimate standard in the absence of exact results for the 
finite wing since consideration of equation (16) will show that this 
eqmtion is at least as convergent when expanded to a certain power of 
03 for a finite wing as when expanded to the same power of 03 for a 
two-dimensional wing. In figure 4 the con^ionents of the section moment- 
curve slopes defined in equations (i«D) are plotted against 03 for 
M = 10/9 and Xq = 0.5* The three dashed-line curves appearing in each 

of the four parts of figure It represent approximate results obtained by 
substitut ing the components given in equations (32) into the appropriate 
expressions of equations (1|0) and retaining terms involving up to 
the third, fifth, and seventh powers . The solid-line curves represent 
exact results obtained by using the moment coefficients tabulated in 
reference 3- Curves of the components of the section lift-curve slope 
are not included because they show essentially the same convergence 
tendencies . 


As may be noted in figure k, the value of o3 at which a particular 
approximation departs from the exact theory is essentially the same 
regardless of the moment conponent considered. By e3ctendlng the 
results of reference 1 ^represented by the curves in fig. k labeled 

"expansion to a?,'' except that in the present paper an term has 

been added to the component 2iM2^ to ap and to 03'^^, the value of 03 

at which the approximation and the exact theory depart at M = IO/9 
has been increased from 0.8 to I.3 and 2.0, respectively. Further inves- 
tigation will show that these values of B remain essentially the same 
regardless of the Mach nvtmber involved. The ranges of reduced frequency k 

(^equal to ^ which the various approximations adequately rep- 


resent the exact theory may therefore be ejcpressed as follows : 


For expansion to a?. 
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and for expansion to 




From these expressions it is readily seen that expanding to rather 

than oB? more than doubles the allowable k range. For exanqple, at 
M = ID/9 the limiting value of k is Increased from O.O76, which is 
too s ma ll, to 0.I9 idilch is sufficiently large to include most flutter 
cases that would likely be encountered at this Mach number. 


In order to study the effect of aspect ratio on the forces and 
moments, consideration is given to the spanwise variation of section 
lift and moment coefficients . The curves in figures 5 and 6 serve as 
illustrations . In figure 5 "the magnitudes of the slopes obtained by 
siJbstituting equations (29) to (32) into equations (llO) are plotted 
against the spanwise variable 1 - ^ for the conditions M = I.3, 

A = 4 , Xq = 0.i<-13, and tu = 0.5Ml-. The quantity ^ is the previously 
used variable y divided by A/2. The phase angles associated with 
the magnitudes of the slopes shown in figure 5 plotted in figure 6 } 
for exan^ile, the phase angle associated with the magnitude 

= 4k^^L^2 + Lij.^ is = tan"^ The portions of the curves 

in these figures in the range O'^l-^^O.ll- (l/p ^ y ^ A/2) also 
apply to a two-dimensional wing under the conditions listed. The effect 
of aspect ratio may be noted in the tip region 0.4<1-^^1.0 
(0 ^ y < l/p) . 


The phase angles 9^, and 9^ in figure 6 are of particular inter- 
est. In keeping with the fact that a component of the force due to h 
in phase with the velocity or a coirponent of the moment due to a 
in phase with the angular velocity d would be destabilizing, values 
of 9ii between 90° and 270° and values of 9am between 0° and l80° 

would indicate possible dynamic Instability in pure bending and pure 
torsion, respectively. The coefficient however, is always positive 

so that 9jj^ is between -90° and 90°; this agrees with the well-known 

fact that pure bending oscillations are always damped, at least in poten- 
tial flow. The phase angle 9am^ "the other hand, is less than l80° 

in the tTO-dimenslonal region and corresponds to \mdamped or unstable 
conditions; whereas in most of the region affected by the wing tip the 
phase angle 9am is greater than l80°, which corresponds to damped or 
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stable conditions. The main feature of figure 6, therefore, is that 
aspect ratio has a stabiliziog effect with regard to the motion a. 

This result seems to indicate that a flutter analysis which tah.es into 
account the spanwise variation of aerodynamic forces would yield a higher 
flutter-speed coefficient than one based on two-dimensional aerodynamic 
forces. 


Upon substituting equations (37) s^d ( 38 ) into equations (^2), the 
overall effect of aspect ratio on the total lift- and moment-curve slopes 
can be calculated for particular values of the parameters M, ffi, Xq, 

A. Although the effect of aspect ratio A may change considerably 
with small changes in one or more of the other parameters, some insight 
into the overall effect may be gained from calculations in which, together 
with A, one of the parameters ffi, Xq, or M is varied while the other 

two remain fixed. Figures 7 8 show the effect on dCLy^da and dCjjj^da 

of varying A and while keeping M and Xq fixed, and figure 9 
shows the effect on of varying A and xq while keeping M 

and as fixed. The effect on the total lift- and moment-curve slopes 
of varying A pmfl M can be extracted from figures 4 and 5 of refer- 
ence 1. 


In figure 7 the magnitudes aufl in figure 8 the associated phase 
angles of the slopes and dCm/da are plotted against cS for 

M = 1 . 3 , Xq = 0 . 5 , and several values of A. In figure 9 the limiting 
value as S — > 0 of the slope dCm^do. is plotted against Xq for 

M = 1.3 and three values of A. The curves in figures 7, 8, and 9 
labeled A = l/p apply to the least value of A permitted by the fore- 
going analysis in that they represent calculations for the combination 
of aspect ratio and Mach nimiber that causes the Mach line from one wing 
tip to intersect the opposite tip at the trailing edge. 

As illu strated in fig\xre 7(s.), a decrease in aspect ratio for a 
given value of o5 produces a decrease in the magnitude of the total 
lift-curve slope. This trend is in general true for the range of £D 
considered »ud is not affected by a change in Xq or M (see fig. 4 

of ref. 1 ) . In the case of the magnltoade of the total moment-curve 
slope, however, the trend with aspect ratio is greatly dependent on Xq. 

As may be seen from figures 7(^) 9 (see also fig. 5 of ref. l) , the 

magnitude may increase, decrease, or perhaps not vary at all (as is the 
case in fig. 9 at xq = O. 67 ) with a decrease in aspect ratio, depen din g 
on the value of Xq under consideration. This seemingly anomalous 
behavior of the total moment-curve slope is associated with the change 
in center of pressure, for prescribed motions of the wing, with change 
in aspect ratio. Also note in figures 7 9 that the departure from 

the two-dimensional (A = cxj) results is rather large for aspect ratios 
less than abo\it it. 
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Prom figure 9 for A = 00 the following phase angles are obtained: 
0AM = 180° for Xq < 0.5 and 0^ = 0° for Xq > 0.5* The phase 

angle of 90° shown in figure 8(b) for A = 00, oi = 0, and Xq = O.5, 
although associated with a slope of zero magnitude and not really signif- 
icant, is therefore a transitional value. Note in figure 8(b) that 0j^ 

is less than 0° for A less than 6 } this result also indicates that 
aspect ratio is stabilizing. 


APPLICATION OF FORCIS AND MOMENTS TO FLUTTER 
Method of Flutter Analysis 


For the purpose of applying the foregoing results, a flutter 
analysis of the Rayleigh type for uniform cantilever wings is now con- 
sidered. Such an analysis involves the selection of a set of modal func- 
tions to approximate the flutter mode, the formation of the flutter deter- 
minant, and the solution of this determinant for the flutter condition. 
(Althou^ the use of the force and moment coefficients derived herein 
for the undlstorted rectangvilar wing results, in certain cases, in a big 
In^jrovement in accuracy over the use of two-dimensional flow coefficients, 
perhaps a more accxirate but much more cimibersome procedure would be to 
use coefficients derived for a distorted wing.) 

The flutter mode of the uniform cantilever wing is assumed to be 
adequately represented by proper combination of the uncoi:5)led first 
bending and first torsion mode shapes of the wing. The flutter deter- 
minant based on these two modes has the form 



-^ha 




(^ 5 ) 


where the determinant elements are given by (see fig. 10 for coordinate ^) 



(¥te) 

(lA-b) 
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Aah = (M^ + 1M2) <3-^ 


(44c) 



and -wiiere co is the flutter frequency. In the event that two-dimensional 
air forces nrii^ moments are used, the force and moment coefficients appear 
as constants in the integrals and can he factored from under the integral 
signs, nri(^ the determinant elements become 



-\th ~ 



4 




Mn + 



Zo^h 


Aoa = 






> 


(1^5) 


(It may he noted that eqs. (44) and (45) are not restricted to any 
particular Mach number range if and Mj[ (i = 1^ 2, 5^ ^) defined 

in accordance with eq. (28) are taken to apply at either subsonic or 
supersonic speed.) The uncotq>led first bending mode shape Zj^ and the 

first torsion mode shape Z^ needed for the evaluation of the aero- 
dynamic integrals of equations (44) are plotted against 1 - ^ in fig- 
ure 11. The mode-shape integrals of equations (44) and (45) are given 
■by 



= 0.25 
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Zh^ta, 


0.357 



0.50 


The flirtter condition is determined from the nontrivial solution 
of equation (lj-3) • This solution he obtained by various means (see 
ch. XIII of ref. 7) . the present paper equation (ii-3) was solved for 
(0^1 0 ^. For a particular wing and Mach number, for which 

Xq, Xqj and 

was varied mtil the ratio 

tion (31) matched that of the wing. In this manner the values of k 
and c%ya) at flutter and consequently the flutter -speed coefficient 

v/bobt were determined for the wing at the selected value of M. 


the imknown 
values of M, 
frequency k 


are specified, the reduced 
determined from equa- 


Flutter Calculations 


The foregoing method of analysis was used to obtain flutter char- 
acteristics in the Mach number range I0/9 % M ^ 1D/6 for a rectangular 

wing with the following properties: A = 4.55^ i/k = 95 *5# Xq = 0.341, 
X(j = 0.350^ = 0.39) anii <%/<%, = 0.5^. CalcxiLations were made on 


the basis of the con^ranents of the section force and moment coefficients 
given in eqmtions (29) to (32). A numerical method for evaluating the 

aerodynamic integrals resulting from the use of these components (for 
exaiaplej J' in eq. (44a) ^ is outlined in appendix B. 


For con5)arison, calculations were also made by using force and moment 
coefficients for a two-dimensional wing (strip theory) . In this latter 
case the determinant elements are of the form given in eqmtlons (45) 
wherein the components and M^ are obtained from either reference 3 

or equations (3I) and (32) of the present paper. In figure 12 the two 
curves resulting from these calculations are plotted in the fonn of 
flutter-speed coefficient "V^bo^ against Mach number M. As may be 

noted in the figure, the cuorves are well-separated at M = I0/9 bxit 
tend to converge as the Mach number is increased. At M = I0/9 the 
flutter speed obtained by using coefficients for a two-dimensional wing 
is about 62 percent of that obtained by using coefficients for a rectan- 
gular wing, whereas at M = 10/6 it is about 95 percent. 
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CONCLUDING ElMAEKB 


The lineeLTized theory for compressible unsteady flow has been xised 
to derive expressions to the seventh power of the frequency of oscilla- 
tion for the velocity potential and section and total forces and moments 
for a thin rectangvilar wing harmonical l y oscillating in pitch and vertical 
translation and moving at supersonic speed. These expressions are exten- 
sions, althoti^ obtained by a different method, of similar expressions 
derived to the third power of the frequency of oscillation in NACA 
Report 1028. As a result of this extension the largest value of reduced 
frequency for which the theory is sufficiently exact has been increased 

from k = to k = ^ where M is the Mach ntimber. 

• m2 m2 

For example, at M = I£)/9 expansion to the third power of the frequency 
is accurate up to aboixfc k = O.O8, whereas expansion to the seventh power 
is accurate up to k = 0.2. The first of these values for k is too 
small but the second is probably large enough to include most rectangular- 
wing flutter cases likely to occur at M = 10/9- 

Flutter calculations were made in the Mach number range from 10 /^ 
to 10/6 for a rectangular wing of aspect ratio 4.53 4y using the section 
force and moment coefficients for a rectangular wing derived in the pres- 
ent paper and the force and moment coefficients for a two-dimensional 
wing given in NACA Report 846. Comparison of the two curves of flutter 
speed indicated, as may be expected, that the rose of finite-wing coef- 
ficients is very influential at Mach numbers near unity but becomes less 
so as the Mach nimber is increased. At M = I0/9 for the particular 
wing analyzed, the flutter speed obtained by using the section coeffi- 
cients for a two-dimensional wing is about 62 percent of that obtained 
by Tosing coefficients for a rectangular wing; whereas at M = IO/6 it 
is aboiit 95 percent. Since fl\jtter speeds calculated on the basis of 
two-dimensional-wing coefficients are generally very conservative in 
the low supersonic speed range, it -vrould. seem that flutter speeds obtained 
by using the section coefficients for rectangular wings of the present 
paper may conpare more favorably with experiment. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., Jantiary l4, 1954. 



lA 


MCA OIN 5076 


55 


APPENDIX A 

SOME IKTEGRAIEED VALUES OP AND 

From eguation (19) the values of the functions F^^ needed in equa- 
tion (22) are as follows : 

“ Vpy<x - 33 t) + X l^ftr/x 

Fg - l/ftr<x - (U) + ^ 

3 1575 5 

P|5 - 6?x? + Topyx^ + iapf^ jaepy*x + y^x - ») + 

p^ . e?i^.->--2?3te5 * gaipapext + ?S>ffV>x2 t yi atVix 4- ypy{x - By) + ^ 

Pg ■ "igxT t t&aAay ■>■ •*• * aoi>fiBV ^ j| 

Pj . a»x8 ■!■ gaftx Tftr + ■roax^P^a <• aj6lBgs35,3 t e9&tA^^+ loaliox?^??? ■!■ : 22 eatg? 63.6 4- ifi^gasTp? ■(. »7fiaB6;,8 ^x - f j) + ^ oto-l\/^ 

The fmctions F^^ to F^ were given previously in reference 1 but are 
repeated here for the sake of coiiQ)leteneBs . 

The values of F^^^ needed in equations (29) and (50) are obtained 

by substituting the value of x at the trailing edge (x = l) into the 
preceding values of F^. That is. 


^n = ^n(l^y) 
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APPENDIX B 

NIMERICAL METHOD OP EVALUATING AERODYNAMIG 
INTEGEIAI5 IN PLUTTER DETERMINANT 


A n\mierical method of Integration for evaluating the aerodynamic 
integrals of eqmtions (1A-) is now outlined. The outline is based on 
the two cases discussed after equation (28) and depicted in figure 3* 

Case 1 (see fig. 3 (a-))j The spanwise stations at which values for 
the components of the section force and moment coefficients given in 
eqtiations (29) to (32) are to he found are shown in figure 10. Since 
equations (29) and (30) reduce, respectively, to equations (3I) and (32) 

for 5 > 2/AP, the values of needed to obtain the coniponents of the 

section force and moment coefficients at the various stations are as 
follows ; 


Station 


y 


1 

h 

?2 

^3 

h 

^5 

h 


?8 

^9 

- Py) 

0 

0 

0 

iAp 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


'IFa 

a/6 

.95661 

.37005 

.217® 

.15332 

.11^ 

.09640 

.08137 

.07040 

.06204 

2 


1/2P 

1/2 

«A 

1.28540 

.55937 

.33958 

.23921 

.1^5 

.14848 

.12464 

.10742 

.09440 

3 

3Xp 

3AP 

VFa 

«/3 

1.48021 

.70402 

.45010 

.32598 

.25329 

.20601 

.17304 

.14562 

.13046 

5, 6 


kL/p 

0 

«/2 

I.570S0 

•7»540 

•52360 

.39270 

.31416 

.26180 

.22440 

.19635 

.17'*53 


Note that these values of do not vary with A and p so long as 

Ap > 2. From the expressions for and in equation (29) ^ the 

preceding table of values of Fn, the curve for in figure 11, and 

the flutter parameters 

k = 0.10 A = 4 

P = 0.8307 H = ^ = 0-1505 


Xo = 0.431 


A = 0.5 - — = 0.1990 

AP 
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the following sample integration table based on Simpson's parabolic 
integrating rnile may be devised: 


Station, n 


Zh^ 

Ij^ + dLTi^ 

In = (Li + 

0 

0 

1.0000 

0 

0 

1 

.1505 

.6275 

.14.549 + 7 .l 684 i 

.2854 + 4.49651 

2 

.3010 

.3469 

•9335 + 9-59321 

.3238 + 3-52811 

5 

.4515 

.1584 

1.3839 + 10.99591 

.2192 + 1.74151 

4 

.6020 

.0520 

1.6646 + 11.62561 

.0865 + . 6 o 42 i 

5 

.8010 

.0041 

1.6646 + 11.62361 

.0068 + .04761 

6 

1.0000 

0 

1 

1.6646 + 11.62361 

0 

/ (Li + IL2) + lUi + 2I2 + 4 lj + I4) + |(ll, + 4I5 + Ig) 

= 0.1456 + 1.66851 


where denotes the value of the integrand at station n» In a 

similar manner the remaining aerodynamic integrals of eq^uations ( 44 ) can 
be evaluated. 

Case 2 (see fig. 5 (b)): In this case, as illustrated by equation ( 55 ) ^ 

the con^jonents of the force and moment coefflclents_for the wing sections 
passing throiigh region contain the functions Fn(A - y) in addi- 

tion to the functions Fn(y) . As a res.ult, there is no arr a ngement of 
spanwise stations which will yield values of that do not vary with 

A and p. For this reason, it seems advisable simply to choose stations 
at even increments along the span for the numerical integration. 
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(b) Case 2. -Mach lines from tips intersect on wing 
but not ahead ot midchord (IIA/3^2). 


Figure 3.- Sketch "based on nondimensional coordinates x and y 
illustrating different Mach line locations accounted for in 
analysis . 
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(a) Real part of moment-curve slope associated with vertical 
translation of wing . 


Expansion to 6)^ 

Expansion to fi)® 

Expansion to cD^ 

Exact (ref. 3 ) 


(b) Imaginary part of moment-curve slope associated with vertical 
translation of wing. 


Figure 4 .- Con5)arison of two-dimensional moment-curve slopes "based on 
exact and approximate theory as a function of frequency parameter S 
for M = 10/9 and Xq = 0.5. 





(d ) Imaginary part of moment-curve slope associated with 
pitching of wing . 


Figure It.- Concluded. 
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Figure Magnitudes of secticm. lift-curve slope and section moment -curve 
slope for rectangular wing of aspect ratio L plotted against spanwise 
position in tenths of semi span meas-ured from wing center for M = 1.5, 
Xq 0.415, and = 0.544. 
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(c) Moment-curve slope associated with vertical 
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Pigure 5-- 
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(d ) Moment-curve slope associated with pitching of wing 
Conclioded. 



(a) Phase angle between lift vector due to vertical 
translation and vertical-velocity vector V\, 
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(b) Phase angle between lift vector due to pitching 
and angular- displacement vector a. 


Figure 6.- Phase angles fcrr rectangular iring of aspect ratio plotted 
againat spanvise position, in tenths of semispan measured froci wing 
center for M = 1.5, Xq = 0.415, ajid S = 0.544. 
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Figure 6.- Concliiied 
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Figure 9-- Total moment-curve slope associated with pitching as a function 
of Xq for M = 1-3, — > 0, and three values of A. 
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Figure 11.- Uncoupled first tending mode slaape and first torsion 

mode shape for a uniform cantilever wing. 
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Figure 12.- Flutter-speed coefficients plotted against Mach number calcu- 
lated from rectangular-wing aerodynamic coefficients and from two- 
dimensional aerodynamic coefficients for A = ^.53, 1/k = 95*3^ 

Xq = 0 . 3 ^ 1 , = 0 . 350 , ro,2 = 0 . 39 , and qt/o^ = O. 583 . 
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